Fluctuations of energy flux in wave turbulence 
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We report that the power driving gravity and capillary wave turbulence in a statistically stationary 
regime displays fluctuations much stronger than its mean value. We show that its probability 
density function (PDF) has a most probable value close to zero and involves two asymmetric roughly 
exponential tails. We understand the qualitative features of the PDF using a simple Langevin type 
model. 
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When a dissipative system is driven in a statistically 
stationary regime by an external forcing, a given amount 
of power per unit mass, e, is transfered from the driv- 
ing device to the system and is ultimately dissipated. In 
fully developed turbulence, a flow is driven at large scales 
and nonlinear interactions transfer kinetic energy toward 
small scales where viscous dissipation takes place. In the 
intermediate range of scales (the inertial range) the key 
role of the energy flux e has been first understood by Kol- 
mogorov [l[. Using dimensional arguments, he derived 
the law E(k) oc e 2 / 3 fc~ 5 / 3 for the energy density E(k) as 
a function of the wavenumber k. Kolmogorov type spec- 
tra have been derived analytically in wave turbulence, 
i.e. in various systems involving an ensemble of weakly 
interacting nonlinear waves (see for instance [2| for a re- 
view). In all cases, it has been assumed that e is a given 
constant parameter. However, it should be kept in mind 
that e is not an input parameter in most experiments or 
simulations of dissipative systems. Its value is not ex- 
ternally controlled but determined by the impedance of 
the system. In addition, as we have already shown for a 
variety of different dissipative systems [H, 0, H| , the en- 
ergy flux or related global quantities, strongly fluctuate 
in time although being averaged in space on the whole 
system or on its boundaries. These fluctuations should 
not be confused with small scale intermittency which oc- 
curs in fully developed turbulence. The later is related to 
the spotness of dissipation in space 0] and its description 
does not involve a time dependent e. 

Here we study the fluctuations of the injected power in 
wave turbulence. Gravity-capillary waves are generated 
on a fluid layer by low frequency random vibrations of 
a wave maker. By measuring the applied force on the 
wave maker and its velocity, we determine the instanta- 
neous power I(t) injected into the fluid. We observe that 
it strongly fluctuates. Its most probable value is 0. rms 
fluctuations oj up to several times the mean value (/) are 
observed, and the probability density function (PDF) of 
/ displays roughly exponential tails for both positive and 
negative values of /. These negative values correspond to 
events for which the random wave field gives back energy 



to the driving device. We show how fluctuations of the in- 
jected power depend on the system size and on the mean 
dissipation and we study their statistical properties. 




FIG. 1: (color online) Time recordings of the velocity of the 
wave maker and the force applied to the wave maker by the 
vibration exciter. The fluid is mercury, with h = 23 mm. 
Both PDF are Gaussian (dashed lines) with zero mean value. 

The experimental setup, described in Q, consists of a 
rectangular plastic vessel, with lateral dimensions 57 x 50 
or 20 x 20 cm 2 , filled with water or mercury (density 13.6 
times larger than water) up to a height, h = 1.8 or 2.3 cm. 
Surface waves are generated by the horizontal motion of 
a rectangular [L x H cm 2 ) plunging PMMA wave maker 
driven by an electromagnetic vibration exciter. We take 
1.2 < L < 25 cm and H — 3.5 cm. The wave maker is 
driven with random noise excitation below 4 or 6 Hz. 

The power injected into the wave field by the wave 
maker is determined as follows. The velocity V(t) of the 
wave maker is measured using a coil placed on the top of 
the vibration exciter. The voltage induced by the mov- 
ing permanent magnet of the vibration exciter is propor- 
tional to V(t). The force Fji(t) applied by the vibration 
exciter on the wave maker is measured by a piezoresistive 
force transducer (FGP 10 daN). The time recordings of 
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V(t) and FA(t) together with their PDFs are displayed 
in Fig. 1. Both V(t) and FA(t) are Gaussian with zero 
mean value. For a given excitation bandwidth, the rms 
value ay of the velocity fluctuations of the wave maker 
is proportional to the driving voltage U applied to the 
electromagnetic shaker and does not depend on the fluid 
density p. On the contrary, the standard deviation ap A 
of the force applied to the wave maker is decreased by the 
density ratio (~ 13) when mercury is replaced by water. 
We have checked that op A oc pSpay where Sp = Lh is 
the immersed area of the wave maker. This linear behav- 
ior has been measured on one decade up to ap A ~ 2 N 
and ay ~ 0.1 m/s. 




FIG. 2: (color online) PDF of I(t)/{I) for mercury: container 
size 57 x 50 cm (red) and 20 x 20 cm 2 (black) (h = 18 mm). 
Dashed lines are the related predictions from equation (f4]) 
without fitting parameter. Inset: time recording of I(t). 

When the wave maker inertia is negligible, the power 
lit) injected into the fluid is roughly given by F>i(t)V(i) 
(see below). The time recording of I it) is shown in the 
inset of Fig. 2. Contrary to the velocity or the force, 
the injected power consists of strong intermittent bursts. 
Although the forcing is statistically stationary, there are 
quiescent periods with a small amount of injected power 
interrupted by bursts where lit) can take both positive 
and negative values. The PDFs of //(/) are displayed 
in Fig. 2. They show that the most probable value of 
I is zero and display two asymmetric exponential tails 
(or stretched exponential in the smaller container). We 
observe that events with I it) < 0, i.e. for which the 
wave field gives back energy to the wave maker, occur 
with a fairly high probability. The standard deviation a/ 
of the injected power is much larger than its mean value 
(7) and rare events with amplitude up to la are also 
detected. Typical values obtained when ay ~ 0.05 m/s 
are ap A ~ 1 N, (I) ~ 30 mW, oj ~ 100 mW for mercury. 
Our measurements also show that aj oc (I) — cpSpa v , 
where c has the dimension of a velocity (c ~ 0.5 m/s and 



slightly increases when the container size is increased) . 

We also observe in Fig. 2 that the probability of neg- 
ative events strongly decreases when the container size 
is increased whereas the positive fluctuations are less af- 
fected. This shows that the backscattering of the energy 
flux from the wave field to the driving device is related to 
the waves reflected by the boundary that can, from time 
to time, drive the wave maker in phase with its motion. 
We note that we have less statistics for the negative tail 
of the PDF when the size of the container is increased. 

We recall that the statistical properties of the fluc- 
tuations of the surface height have been studied in 0]: 
they involve a large distribution of amplitude fluctua- 
tions. Their frequency spectrum is broad band and can 
be fitted by two power laws in the gravity and capillary 
regimes. The power law exponent in the capillary range 
is in agreement with theoretical predictions. The one in 
the gravity range depends on the forcing, as also shown 
in [8( . The scaling of the spectrum with respect to the 
mean energy flux (I) is different from the theoretical pre- 
diction both in the gravity and capillary ranges. These 
discrepancies can be ascribed to finite size effects 0, IH ■ 
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FIG. 3: (color online) Effect of the inertia of the wave maker: 
PDF of F A V (black) and of (F A — M V)V (red) for mercury 
(h = 23 mm) (left) and for water (h = 23 mm) (right). Using 
F A V to estimate I leads to an error on the standard deviation 
oi that is less than 5% for mercury but that reaches 50% for 
water. 

We first emphasize the bias that can result from the 
system inertia when one tries a direct measurement of the 
fluctuations of injected power. The equation of motion 
of the wave maker is 

MV = F A (t) + F R (t), (1) 

where M is the mass of the wave maker and Fp(t) is 
the force due to the fluid motion (V = dV/dt). The 
power injected into the fluid by the wave maker is I(t) = 
-F R (t)V(t). When MV is not negligible, I(t) gener- 
ally differs from FA(t)V(t) which is experimentally deter- 
mined. This obviously does not affect the mean value (I) 
but may lead to wrong estimates of fluctuations. Using 
an accelerometer, we have checked that M V is negligi- 
ble compared to Fa when the working fluid is mercury. 
This is shown is Fig. 3 (left) where the PDF of FaV and 
FpV = (Fa — MV)V are compared. On the contrary, 
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inertia cannot be neglected for experiments in water for 
which an error as large as one order of magnitude can be 
made on the probability of rare events if one use F^V 
to estimate I (right). Thus, the correction due to MV 
has been taken into account in water experiments. There 
exist only a few previous direct measurements of injected 
power in turbulent flows and those type of inertial bias 
have never been taken into account 



Fluid effect 



I/M,F) and integrating over F gives 




FIG. 4: (color online) Effect of fluid properties on the PDF 
of the injected power: (red) mercury; (black): water (h — 18 
mm; 20x20 cm 2 container). Solid lines indicate the value of 
(I) 

The PDFs of injected power for the same driving in 
the same container for water and mercury are displayed 
in Fig. 4. The asymmetry of the PDF is much larger 
with mercury. This is related to its larger mean energy 
flux, i.e., mean dissipation, as shown below. 

The qualitative features of the PDF of injected power 
can be described with the following simple model. 
Guided by our experimental observation of the linearity 
of a f a m a v > we assume that the force Fr due to the fluid 
can be roughly approximated by a friction force —M"fV 
where 7 is a constant (the inverse of the damping time of 
the wave maker) . We are aware that a better approxima- 
tion to the force due to the fluid should involve both V 
and an integral of V(t') with an appropriate kernel. Thus 
we only claim here to give a heuristic understanding of 
the qualitative properties of the PDF of /. Modelling the 
forcing with an Ornstein-Uhlenbeck process, we obtain 



V = -7V + F, F = -pF + £, 



(2) 



where (3 is the inverse of the correlation time of the ap- 
plied force (F = Fa/M) and £ is a Gaussian white noise 
with (£(t)£(f)) = DS(t - t'). The PDF P(V,F) is the 
bivariate normal distribution T3, U | 
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where Kq(X) is the zeroth order modified Bessel func- 
tion of the second kind. Using the method of steep- 
est descent, this predicts exponential tails, P(X) ~ 



(l/y/\X\)exp(rX-\X\) where X = I/[(l-r 2 )a v a F ]. In 
addition, we have (I) = D/[2f3(j + f3)\ — rayap- Thus, 
(j4]) is determined once (I), ay and of have been mea- 
sured and can be compared to the experimental PDF 
without using any fitting parameter. This is displayed 
with dashed lines in Fig. 2. Taking into account the 
strong approximation made in the above model, we ob- 
serve a good agreement in the larger container. More 
importantly, this model captures the qualitative features 
of the PDF: its maximum for I — and the asym- 
metry of the tails that is governed by the parameter 
r = VlJW+Wj = (I)/{owf a )- For given a v and <t Fa , 
the larger is the mean energy flux, i.e., the dissipation, 
the more asymmetric is the PDF. For mercury, direct de- 
termination of r from the measurement of (I), <rv and 
c f a gives r ~ 0.7 for the large container and r ~ 0.6 for 
the small one, in qualitative agreement with the differ- 
ent asymmetry of the PDF in Fig. 2. Smaller values of r 
are achieved in water for which the dissipation is smaller. 
The PDFs are more stretched for water in particular in 
the smaller container. 
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FIG. 5: (color online) PDFs of the injected power I T averaged 
on a time interval r: r = 1, 3, 11 and 50r c , where r c = 0.03 s 
is the correlation time of I(t). Solid lines indicate the value 
of (J) (water, h = 23 mm). 

We now consider the injected power averaged on a time 
interval r 

Ir(t) = ~ I(t')dt'. (5) 
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The PDFs of I T for t/t c = 1, 3, 11 and 50 where r c is the 
correlation time of I(t), are displayed in Fig. 5. They 
become more and more peaked around I T ~ (I) as they 
should. However, one needs to average on a rather large 
time interval (r ~ 50t c ) in order to get a maximum prob- 
ability P(I T ) for I T = (I) (Fig. 5, bottom right). Then, 
the probability of negative events become so small that 
almost none can be observed. Fig. 6 shows that the 
quantity, i log p7~y~77TO f° r different values of t that 
has been predicted to be linear in I T /(I) when the hy- 
pothesis of the fluctuation theorem (in particular time 
reversibility) are fulfilled [HI, [l3| ■ As we clearly observe 
in Fig. 6, this is not the case in general for dissipative 
systems. As already mentioned Q and studied in details 
14j . the linear behavior reported in several experiments 



or numerical simulations results from the too small values 
of I T /(I) that are probed when r ^> r c . Large enough 
values are obtained in the present experiment and the 
expected nonlinear behavior is thus reached. The shape 
of the curve in Fig. 6 is found in good agreement with 
the analytical calculation [loT ] performed with a Langevin 
type equation with white noise. 
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FIG. 6: (color online) Plot of i log ^j/)^ for 16 <t/t c < 
39 [t/t c = 17 (*), 19.5 (o), 22 (□), 25 "(<>), 28 (pentagram), 
30.5 (V), 33.5 (hexagram), 39 (<)]. Langevin model of ref- 
erence [15[ : 47 for e — I T /{I) < 1/3 (dashed line) and 
77e/4 + 37/2 - 7/(4e) for e > 1/3 (solid line) with 7 = 5 
Hz. 



Finally, we emphasize that a fluctuating injected power 
implies fluctuations of the energy flux at all wave num- 
bers in the energy cascade from injection to dissipation. 
In any system where an energy flux cascades from the in- 
jected power at large scales to dissipation at small scales, 
one has for the energy _E< for wave numbers smaller than 
k within the inertial range, E < = I(t) — &(k,t) = R, 
where $(fc,i) is the energy flux at k toward large wave 
numbers. Thus ' (R(t)R(0))<1t = in order to prevent 
the divergence of (E^). Dimensionaly, this implies that 



a%Tk does not depend on k where cr$ is the standard 
deviation of the energy flux and Tk is its correlation time. 
If this dimensional scaling is correct, fluctuations of the 
energy flux are expected to increase during the cascade 
from large to small scales since Tk decreases (for instance, 
cr$ oc fc 1 / 3 for hydrodynamic turbulence). Such fluctua- 
tions have been found numerically and experimentally 
in hydrodynamic turbulence [lij ]. To which extent, this 
is related or modified by small scale intermittency [13] 
remains an open question. 
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